Comparative study for non-statistical fluctuation of 
net- proton, baryon, and charge multiplicities 
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We calculate the real and non-statistical higher moment excitation functions (-^/sjvjv = 11.5 to 
200 GeV) for the net-proton, net-baryon, and the net-charge number event distributions in the 
relativistic Au+Au collisions with the parton and hadron cascade model PACIAE. It turned out 
that because of the statistical fluctuation dominance it is very hard to see signature of the CP 
singularity in the real higher moment excitation functions. It is found that the property of higher 
moment excitation functions are significantly dependent on the window size, and hence the CP 
signatures may show only in a definite window for a given conserved observable non-statistical 
higher moments. But for a given widow size, the CP singularity may show only in the non-statistical 
higher moment excitation functions of a definite conserved observable. 
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I. INTRODUCTION 

One fundamental aim of relativistic heavy-ion colli- 
sions is to explore the phase transition from the hadronic 
matter (HM) to the quark-gluon matter (QGM). Ev- 
idences for the strongly coupled quark-gluon plasma 
(sQGP) were reported years ago [Irtl]. The RHIC beam 
energy scan (BES) was then proposed [f|. Attempts arc 
being made to experimentally locate the QCD critical 
point (CP), where the first order phase transition gives 
over to the "crossover" Q in the QCD phase diagram of 
temperature T vs. baryon chemical potential hb @~G3 ■ 
Because a cms energy (^snn) in the heavy-ion collisions 
maps out & fiB value, a T value can be extracted by fitting 
the particle ratio in the hadronic resonance gas (HRG) 
model [Ulllli to corresponding data in the heavy-ion col- 
lisions. The scanning over the phase diagram can then 
be performed. Meanwhile, the search for CP signatures 
could also be done at each cms energy via the singularity 
(nonmonotonicity) shown in higher moment excitation 
functions (higher moment as a function of reaction en- 
ergy) of the conserved observables. Since the birth of 
RHIC BES considerable progresses have been made in 
this field both experimentally [l3l - fl(| and theoretically 
p7l - l23 |. A upsurge was reached last year at the 7 th In- 
ternational Workshop on Critical Point and Onset of De- 
confinement, Nov. 7-11 Wuhan, China. Although many 
relevant fundamental physics have been demonstrated, 
the mysteries remain. 

Early in 2009 one employed the UrQMD model in- 



vestigating the net-proton, net-baryon, and net-charge 
kurtosis in the central Pb+Pb/ Au+Au collisions from 
E lab =2A GeV to ^5^=200 GeV [H. They did not 
find any evidence of CP singularity in the kurtosis exci- 
tation functions. However, they did find that the kurto- 
sis excitation functions behave quite different among the 
net-proton, net-baryon, and net-charge: The net-proton 
kurtosis becomes slightly negative at low ^/snn and the 
net-baryon kurtosis decreases to a large negative value 
while the net-charge kurtosis keeps zero approximately. 
Similarly, the kurtosis as a function of the rapidity win- 
dow size were also different: Both the net-proton kurtosis 
and net-charge kurtosis are differently fluctuated around 
zero, but the net-baryon kurtosis reaches a large negative 
value at the lager window. On the contrary, the opinion 
in [26| is quite opposite. By assuming that the num- 
ber distribution of a conserved observable is a binomial 
function and the strange baryon can be approximately 
considered, they analytically derived a relation between 
baryon number cumulant and the proton number cumu- 
lant. This means that the net-baryon cumulant has a 
similar behavior to the net-proton. 



*zhoudm@phy.ccnu.edu.cn/zdm@iopp. ccnu.edu.cn 



Recently it has been reported in [2JJ that the effect of 
finite acceptance dramatically influences the cumulants 
of conserved observables, which makes the interpretation 
of experimental data measured in a given acceptance very 
uncertain. In 27] the concept of "required acceptance" is 
introduced, where all the particles, which are necessary 
for the physics in question, are covered. The required ac- 
ceptance is not the full acceptance where all particles in 
the final hadronic state are covered. The authors assume 
that the bias of a real acceptance in a real experiment 
from the required acceptance can be represented by a 
binomial probability distribution, then they derive the 
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general formula that relate the true cumulants which re- 
flect the full dynamics of the system (require acceptance) 
with the actually measured cumulants (real acceptance). 
This relation involves additional moments which cannot 
be expressed by the cumulants and should be simulta- 
neously measured in order to extract information about 
CP singularity. But how to measure these additional mo- 
ments is not yet known. 

We attempt, in this paper, to further explore the rel- 
evant issues using the parton and hadron cascade model 
PACIAE [28} . The net-proton, net-baryon, and the net- 
charge number event distributions and their higher mo- 
ment excitation functions (*/sjvw=11.5 to 200 GeV) in 
the 0-5% most central Au+Au collision are calculated in 
different transverse momentum and pseudorapidity win- 
dows by the PACIAE model. As pointed out in [29| that 
the study of higher moment excitation function singu- 
larity is a matter of dynamical fluctuation, but the sta- 
tistical fluctuation is always dominant. In order to pro- 
nounce the effect of dynamical fluctuation, non-statistical 
higher moments (cumulants) of conserved observables are 
proposed in [2{|. Therefore, we calculate simultaneously 
both the real and the non-statistical higher moment ex- 
citation functions, and comparatively analyze the depen- 
dence of these excitation functions on the window size. 



II. MODELS 

Suppose that the event distribution of the number of 
conserved observable x is P(x), then its n*' 1 moment 
about zero is 



(x {n) ) = J x n P{x)dx, (1) 

and its n th moment about the mean (M = (a;' 1 )) = (x)) 
is 

Af = ((a- - (a;))") = [ (x- (x)) n P(x)dx. (2) 



Therefore its higher moments and moment products in- 
vestigated widely are: 



events are generated by combining particles selected ran- 
domly from different real events, while reproducing the 
event multiplicity distribution of the real events [30]. If 
the ii th real (mixed) moment about the mean calculated 
with real (mixed) events is denoted by (M^p), then 
the corresponding n th non-statistical moment is 
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(8) 

like in (30|. 

PACIAE is based on PYTHIA [U which is devised 
for the high energy hadron-hadron (hh) collisions. In the 
PYTHIA model, a hh collision is first decomposed into 
the parton-parton collisions. The hard parton-parton in- 
teraction is dealt with lowest leading order perturbative 
QCD (LO-pQCD) and the soft parton-parton interac- 
tion is considered empirically. Because the consideration 
of the initial- and final- state parton showers as well as 
the parton multiple interactions, the consequence of a 
hh collision is a partonic state composed of quarks (anti- 
quarks), diquarks (anti-diquarks) and gluons. It is fol- 
lowed by the string construction and fragmentation. So 
a final hadronic state is obtained for a hh collision even- 
tually. 

In the PACIAE model ,28] , a nucleus- nucleus collision 
is first decomposed into nucleon-nucleon (NN) collisions 
according to the collision geometry and NN total cross 
section. Each NN collision is dealt by PYTHIA with 
string fragmentation switched-off and the diquarks (anti- 
diquarks) broken into quark pairs (anti-quark pairs). 
When all NN collision pairs are exhausted, a partonic 
initial state (composed of quarks, antiquarks, and glu- 
ons) is obtained for a nucleus-nucleus collision. This 
partonic initial stage is followed by a parton evolution 
stage, where parton rescattering is performed by the 
Monte Carlo method with 2—^2 LO-pQCD cross sec- 
tions [32]. The hadronization stage follows the parton 
rescattering stage. The Lund string fragmentation model 
and a phenomcnological coalescence model are provided 
for the hadronization. Then the rescattering among the 
produced hadrons is dealt with the usual two body col- 
lision model (28|. In this hadronic evolution stage, only 
the rescatterings among 7T, K, p, n, p(u>), A, A, E, S, SI, 
and their antiparticles are considered for simplicity. We 
refer to [28j for the details. 



skewness: S = M (3) /(Af (2) ) 3/2 , (4) 

kurtosis: k = M< 4 ty(M< a >) a - 3, (5) 

Sa = M (3) /M (2) , (6) 

and 

Ka 2 =MW/MM-3MW. (7) 

We employ the parton and hadron cascade model PA- 
CIAE to generate the real event distributions. The mixed 



III. RESULTS 

We use PACIAE to generate the net-proton, net- 
baryon and the net-charge number event distributions 
in 0-5% most central Au+Au collisions at v / stvw=11.5 
to 200 GeV for the different acceptance sizes. The real 
and non-statistical higher moment excitation functions 
are then calculated, with all model parameters fixed. 
Baryons are identified with p + n + A and their anti- 
particles (the non-strange baryons heavier than neutron 
and the strange baryons heavier than A are assumed de- 
cay already) and the charge with , K ± , p and p. 
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FIG. 1: (color online) Real higher moment excitation func- 
tions of the net-proton (top), net-baryon (middle), and the 
net-charge (bottom) in the 0-5% most central Au+Au colli- 
sions calculated by PACIAE for different window sizes. 



We summarize our results as follows: 

Shown in the Figure [TJ from the top to the bottom, 
are respectively the net-proton, net-baryon, and the net- 
charge real higher moment excitation functions, with the 
mean M, the variance a, skewness S 7 kurtosis re, and the 
products Sa and red 2 displayed in the panels (a), (b), 
(c), (d), (e) and (f), respectively. In the figure, the solid 
circles are the results in the acceptances of 0.4 < px < 0.8 
GeV/c and \y\ < 0.5, the open circles are the results in 
0.4 < pt < 0.8 and \y\ < 1, and the diamonds are the 



FIG. 2: (color online) Real higher moment excitation func- 
tions of the net-proton, net-baryon, and the net-charge in the 
0-5% most central Au+Au collisions calculated by PACIAE 
for a given window size of 0.2 < pr < 1.6 GeV/c and \y\ < 1. 



results in 0.2 < px < 1.6 and \y\ < 1. 

Figure [2] gives, respectively, the net-proton (circles), 
net-baryon (squares), and the net-charge (triangles) real 
higher moment excitation functions for the mean (panel 
(a)), g ((b)), skewness ((c)), kurtosis ((d)), Sa ((e)), and 
K(T 2 ((f)) calculated in the 0.2 < px < 1.6 GeV/c and 
\y\ < 1 acceptances. 

We see from these figures that: 

• The shape of net-proton (net-baryon, net-charge) 
real higher moment excitation function changing 
with the window sizes is approximately similar only 
for the mean and a, see top panel (a) and (b) (mid- 
dle panel (a) and (b), as well as bottom panel (a) 
and (b)) of Fig. [1] This change becomes rather 
complicated for the skewness S, kurtosis re, and the 
product Sa and recr 2 . 

• We see in Fig. [1] that the net-proton kurtosis is 
always positive (cf. top panel (d)) unlike in [25j], the 
net-baryon kurtosis decreases to a small negative 
value at low ^/s^n (cf. middle panel (d)) but not 
to a lager negative value as mentioned in [251 ) , while 
the net-charge kurtosis approximately keeps zero 
(cf. bottom panel (d)) like in [25j ■ 

• It is hard to see any signature, which may indicate 
the CP singularity, from the real higher moment 
excitation functions shown in these figures. Only 
the net-charge real excitation function of kurtosis in 
the window of 0.4 < p T < 0.8 GeV/c and \y\ < 0.5 
as well as 0.4 < px < 0.8 GeV/c and \y\ < 1 may 
indicate some signature (cf. bottom panel (d) in 
Fig. UK. 

Figure [3] gives, from the top to the bottom, the net- 
proton, net-baryon, and the net-charge non-statistical 
higher moment excitation functions for the skewness 
(panel (a)), kurtosis (panel (b)), Sa (panel (c)), and rea 2 
(panel (d)). In the figure the solid circles, open circles, 
and the diamonds are, respectively, the results in the 
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FIG. 4: (color online) Non-statistical higher moment excita- 
tion functions of the net-proton (top), net-baryon (middle), 
and the net-charge (bottom) in the 0-5% most central Au+Au 
collisions calculated by PACIAE for a given window sizes of 
0.2 < p T < 1.6 GeV/c and \y\ < 1. 



energy region may serve as a signature of the CP sin- 
gularity. As for the net-baryon non-statistical excitation 
functions, it is the kurtosis and na 2 which may display 
the CP singularity only in the window of 0.4 < px < 0.8 
GeV/c and \y\ < 0.5, as shown in the middle panels (b) 
and (d). However, one can not find any possible signature 
of the CP singularity in the bottom panels. 

Figure [5] gives, respectively, the net-proton (circles), 
net-baryon (squares), and the net-charge (triangles) non- 
statistical excitation functions for the skewness (panel 
(a)), kurtosis (panel (b)), Sa (panel (c)), and na 2 (panel 
(d)) calculated in the acceptance of 0.2 < p T < 1.6 
GeV/c and |y| < 1. We see in FigH] (a) and (c) that the 
CP singularity may be shown in the non-statistical exci- 
tation functions for the net-proton but not for the net- 
baryon and net-charge. However, a CP signature may 
appear in the non-statistical excitation functions of kur- 
tosis for the net-baryon but not for the net-proton and 
net-charge, as shown in Fig[4](b). 



IV. CONCLUSIONS 



FIG. 3: (color online) Non-statistical higher moment excita- 
tion functions of the net-proton (top), net-baryon (middle), 
and the net-charge (bottom) in the 0-5% most central Au+Au 
collisions calculated by PACIAE for different window sizes. 



acceptances of 0.4 < pr < 0.8 GeV/c and \y\ < 0.5, 
0.4 < p T < 0.8 and \y\ < 1, and 0.2 < p T < 1.6 and 
\y\ < 1. One sees an interesting structure, as shown in 
the net-proton non-statistical excitation functions of the 
skewness and Sa (cf. top panels (a) and (c) in Fig. [3]), 
which becomes more pronounced with increasing the win- 
dow size. The structure that the net-proton skewness and 
Sa go down to a sizably negative value within a small 



We have calculated the real and non-statistical higher 
moment excitation functions (^/s~n~n = H-5 to 200 GeV) 
for the net-proton, net-baryon, and the net-charge num- 
ber event distributions in the relativistic Au+Au colli- 
sions with the parton and hadron cascade model PA- 
CIAE The real event is generated with the PACIAE 
model, while the mixed event is constructed by combin- 
ing particles selected from the real events randomly. The 
non-statistical moment is derived as the difference be- 
tween real moment and the mixed moment. It turned 
out that because of the statistical fluctuation dominance 
29] it is very hard to see any signature of the CP singu- 
larity in the real higher moment excitation functions of 
the net-proton, net-baryon, and the net-charge, which is 
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consistent with the results from the real higher moment 
studies in [Him. 

The work suggests that the non-statistical moment of 
conserved observables may play a potential role in ex- 
ploring the CP singularity in the relativistic heavy ion 
collisions. However, it is found that the property of 
higher moment excitation functions is significantly de- 
pendent on the window size. For a given conserved ob- 
servable, a signature of the CP singularity may appear 
in its non-statistical higher moment excitation functions 
only in a definite window size. But for a given widow 
size, the CP singularity may be shown only in the non- 



statistical higher moment excitation functions of a defi- 
nite conserved observable as the multiplicity distribution 
over rapidity and transverse momentum among different 
conserved observables can be rather different. 
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